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PRIMER  NOTE 

Isolation  of  polymorphic  microsatellite  markers  from  the 
malaria  vector  Anopheles  marajoara  (Diptera:  Culicidae) 

CONG  LI/  RICHARD  C.  WILKERSON  *§  and  DINA  M.  FONSECA  + 

*Department  of  Entomology,  Division  of  Com?mmicable  Diseases  and  Immunology,  Walter  Reed  Army  Institute  of  Research,  503  Robert 
Grant  Avenue,  Silver  Spring,  MD  20910-7500,  t  Academy  of  Natural  Sciences,  1900  Benjamin  Franklin  Parkway,  Philadelphia, 

PA  19103 


Abstract 

Microsatellite-containing  regions  were  isolated  and  characterized  in  Anopheles  ( Nysso - 
rhynehus)  marajoara,  a  primary  vector  of  malaria  parasites  in  northeastern  Amazonia, 
Brazil.  An  enrichment  protocol  yielded  500  positive  clones.  We  designed  primers  to 
amplify  40  unique  microsatellites,  11  of  which  amplified  cleanly  and  were  polymorphic.  A 
survey  of  323  individuals  showed  that  these  loci  are  highly  variable  with  the  number  of 
alleles  ranging  from  11  to  52,  and  expected  heterozygosity  ranging  from  0.64  to  0.95.  These 
markers  will  be  useful  for  studies  of  population  structure  and  intraspecific  variation  in 
A.  marajoara. 
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Anopheles  ( Nyssorhynchus )  marajoara  Galvao  &  Damasceno, 
ranges  from  southern  Sao  Paulo  State,  Brazil  to  Costa 
Rica,  and  was  recently  recognized  as  the  primary  vector 
of  malaria  parasites  in  northeastern  Amazonia,  Brazil 
(Conn  et  al.  2002).  An.  marajoara  is  a  member  of  the  albitarsis 
complex  which  includes  four  largely  isomorphic  species 
that  can  only  reliably  be  separated  using  RAPD  (random 
amplified  polymorphic  DNA)  markers  (Wilkerson  et  al. 
1995a,  1995b)  or  by  PCR  (polymerase  chain  reaction)  of 
diagnostic  rDNA  ITS2  sequence  (Li  et  al.  unpublished). 
Because  of  the  significant  vector  status  of  An.  marajoara, 
its  wide  geographical  range,  and  questions  regarding  its 
taxonomic  status,  we  developed  microsatellite  markers  to 
study  population  and  species  level  questions.  This  is  only 
the  second  Neotropical  malaria  vector  for  which  micro¬ 
satellite  markers  have  been  developed,  the  other  being 
Anopheles  ( Nyssorhynchus )  darlingi  Root  (Conn  et  al.  2001). 

Microsatellite  loci  were  isolated  using  an  enrichment 
protocol  following  Hamilton  etal.  (1999)  with  modifica¬ 
tions  as  in  Keyghobadi  etal.  (2004).  Genomic  DNA  was 
extracted  from  three  specimens  of  An.  marajoara  by  phenol- 
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chloroform  extraction  as  in  Wilkerson  et  al.  (1993)  yielding 
approximately  two  pg  of  DNA  which  was  digested  with 
Haelll,  Msel,  and  Nhel  restriction  enzymes  (New  England 
Biolabs  —  NEB).  The  resulting  200-1000  bp  fragments  were 
blunt-ended,  dephosphorylated,  and  ligated  to  linkers  as 
in  Keyghobadi  et  al.  (2004).  This  DNA  was  then  'enriched' 
for  GT,  GA,  CAC,  GTC,  GGT,  and  GCT  repeats  using 
streptavidin-coated  magnetic  beads  (Dynabeads),  and  cloned 
(Hamilton  et  al.  1999).  Colonies  were  lifted  with  uncharged 
nylon  membranes.  DNA  was  fixed  to  the  membranes,  which 
were  then  hybridized  to  biotin-labelled  (GT)15  (GA)15  (CAC)10 
(GTC)10  (GGT)1C  and  (GCT)10  oligonucleotides  at  65  °C, 
and  positive  colonies  were  detected  using  the  Phototope- 
Star  Chemiluminescent  Detection  Kit  (NEB).  Approximately 
2000  colonies  were  screened,  of  which  about  one  fourth 
were  positive,  and  150  positive  colonies  were  picked,  boiled 
in  100  pL  TE  for  10  min,  and  vortexed.  After  centrifugation 
for  5  min  at  15  000  rpm,  1  pL  of  supernatant  was  PCR 
amplified  with  T3  (5'-AATT A ACCCTCACTAAAGGG-3') 
and  T7  (5'-GTAATACGACTCACTATAGGGC-3')  vector 
primers  and  directly  sequenced  using  Big  Dye  3.0  (Applied 
Biosystems  —  AB).  Fragment  analysis  was  conducted  on  a 
capillary  automatic  sequencer  ABI3100  (AB).  A  blue  colony 
(no  insert)  was  used  as  a  negative  control  in  each  PCR  reaction. 

For  each  microsatellite  sequence  with  sufficient  flanking 
sequence,  flanking  primers  were  designed  with  primer  3 
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Table  1  Characterization  of  11  microsatellite  loci  in  Anopheles  marajoara 


Locus 

Primer  Sequence  (5-3') 

Repeat  Motif 

rot 

No.  of 
allelesf: 

Allele  size 
range  (bp) 

«o 

Gen  Bank 
accession  No. 

M10-4 

F:  HEX-CACirrCCGCTTCACTCCTC’1, 

R:  AACAGTTCGGCGGTCAAATC 

(CA)2, 

56 

52 

211-265 

0.79 

0.95 

AY741201 

Ml -3 

F:  TTTCATAATTATCCCGAAACGAA 

FI*:  CATAATTATCCCGAAACGAACTG 

R:  FAM-GGCAACGATCTGCACTCTTAC 

(ggt)7 

50 

15 

196-221 

0.58 

0.83 

AY741202 

Mll-4 

F:  GTCCTCACGACGACAACGAT 

R:  NED-GTGTTACCAGCGATGGTGTG 

(tga)5 

50 

16 

195-226 

0.51 

0.68 

AY741203 

M5-7 

F:  HEX-CAATTGGTGAACGGTTTCCT 

R:  CGACAGGGTGGTTAGTGGTT 

(tgc)7 

50 

14 

183-215 

0.64 

0.84 

AY741204 

M6-5 

F:  HEX-TGTTTGTCCCACAGCAGAAG 

R:  GCAGCAGATGCCATAAATCA 

(tcg)7  +  (tcg)3 

50 

13 

132-172 

0.69 

0.81 

AY741205 

M2 

F:  ACCAGGGATCAGGGATCGT 

FI:  ACCAGGGATCAGGGATCGTC§ 

R:  FAM-CACCGAGATCCTGCCGTA 

(ggc)5  +  (gac)6 

50 

12 

145-175 

0.29 

0.84 

AY741206 

M7-1 

F:  FAM-ctggaaacctgaaacctgga 

R:  ACGGAACCCTCICTTTGTGA 

(AC)13 

50 

17 

141-169 

0.59 

0.83 

AY741207 

Ml  0-1 

F:  NED-GCTGCCAATGAAAGATTGCT 

R:  GTGTGCGTACACGAGCTACC 

(AC)n 

56 

25 

133-167 

0.67 

0.79 

AY741208 

M8-2-2 

F:  FAM-CAGGAGATGCCACAAATCAG 

R:  AGAGACCAACTTTCGGATGC 

(CGA)a 

45 

13 

105-144 

070 

0.82 

AY741209 

M6-8 

F:  FAM-ATAGACCCCATTCGTCATCG 

R:  TAGGTACGTCCATGGTAGCAGT 

(TGC)10 

50 

11 

100-133 

0.56 

0.64 

AY741210 

M9-2 

F:  CATAGACTTCACACGTACGCACT 

R:  FAM-AATCGGCACTTCCGCTTC 

(GCA)7 

50 

11 

87-119 

0.42 

0.66 

AY741211 

*F  indicates  forward  primer  and  R  indicates  reverse  primer;  +7'a  is  the  touchdown  annealing  temperature  (touchdown  protocol  should  start 
with  a  T  =  T  +  10  °C);  ^Number  of  alleles,  allele  size  range,  and  observed  and  expected  heterozygosity  (H0  and  HE,  respectively),  were 
estimated  from  a  sample  of  323  individuals  from  northern  Brazil;  §altemative  primers  for  populations  with  null  alleles. 


(Rozen  &  Skaletsky  2000).  DNA  of  individual  mosquitoes 
for  population  studies  was  extracted  as  above  (Wilkerson 
et  al.  1993).  The  polymorphism  of  loci  was  tested  by  sizing 
PCR  products  from  12  individuals  (selected  from  different 
localities)  on  6.0%  acrylamide  gels  (Protean  2  system,  BIO¬ 
RAD)  and  visualized  with  ethidium  bromide. 

From  the  40  loci  for  which  primers  were  designed,  we 
selected  1 1  polymorphic  and  easily  scored  loci  for  studying 
population  differentiation  (Table  1).  One  of  the  primers 
on  each  pair  was  labelled  with  FAM,  HEX,  or  NED  (AB). 
Touch-down  PCR  was  conducted  in  a  20  pL  reaction  con¬ 
taining  1  x  buffer  (10  mM  Tris-HCl,  pH  8.3,  50  mM  KC1), 
2.5  mM  MgCl2,  200  |i.\i  dNTP,  5  pmol  each  primer,  0.5  U 
Taq  polymerase  (AB)  and  0.5  pL  DNA  template  (~5  ng). 
PCRs  were  performed  on  a  GeneAmp  PCR  System  9700 
(AB)  with  the  following  conditions:  2  min  at  95  °C  fol¬ 
lowed  by  14  cycles  using  a  touchdown  protocol  (95  °C  for 
30  s,  30  s  at  a  temperature  that  decreased  by  0.7  °C  every 
second  cycle  to  a  'touchdown'  annealing  temperature  (Ta 
in  Table  1),  72  °C  for  30  s),  followed  by  25  cycles  with  a 
standard  protocol  (95  °C  for  30  s,  Ta  for  30  s,  72  °C  for 
30  s),  closing  with  10  min  at  72  °C  to  reduce  stutter  bands. 
Products  from  three  to  six  microsatellite  loci  were  mixed 


and  run  on  an  ABI3100  (AB).  Data  were  automatically 
collected  and  analysed  by  genotyper  software  2.5  (AB). 

Linkage  disequilibrium  for  all  pairs  of  loci  were  tested 
with  genepop  3.3  (Raymond  &  Rousset  1995).  After  Bon- 
ferroni  correction,  one  pair  of  loci  was  found  to  be  signi¬ 
ficantly  linked:  M6-5  with  M8-2-2,  We  compared  observed 
and  expected  frequencies  under  Hardy- Weinberg  equilib¬ 
rium  using  arlequin  (Schneider  et  al.  2000)  with  100  000 
steps  in  the  Markov  chain  and  1000  dememorization  steps. 
After  Bonferroni  correction,  all  loci  conformed  to  Hardy- 
Weinberg  frequencies  except  for  M2.  Upon  further  testing 
we  concluded  this  locus  is  sex  linked  (occurs  on  the  X- 
chromosome)  as  deviations  from  expected  frequencies 
only  occur  in  samples  where  both  males  and  females  were 
examined  (as  opposed  to  only  females). 
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